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We investigate the possibility that the Universe is significantly reionized by the decay products 
of heavy particles. The ionization produced by decay particles implies a high optical depth even if 
the maximum level of ionization ever produced is low (10~^). As a consequence, a high ionization 
fraction {x ~ 0.5) at high redshifts {z ~ 20) fails to fit the cosmic microwave background (CMB) 
spectra at I > 30. Recent CMB data limits the primordial abundance of the decaying particles, 
favoring long decay times. Other significant sources of reionization are still needed at 2; ~ 13. The 
decay process heats up the medium, bringing the expected y distortion to unobservable levels. 

PACS numbers: 95.35.+d, 98.80.Cq 



1 ) Introduction 

The recent results of the WMAP satellite I"] have evi- 
denced that the universe presents a higher optical depth 
that previously anticipated from the observation of ab- 
sorption lines toward quasars at redshift z ~ 6. The 
data suggest that the Universe is highly ionized at red- 
shifts as high as z ~ 20. Since this discovery, scientists 
have investigated different possibilities for producing this 
high-redshift reionization. In the standard reionization 
scenario where the first stars are fully responsible for 
reionization, the WMAP result would imply great difh- 
culties for a number of models, namely the ones in which 
the matter power spectrum at small scale is suppressed. 
Among these are, for example, the warm dark matter 
(WDM) model of structure formation 0, 0] and the cold 
dark matter model with a running spectral index. These 
models wouldn't provide the sufficient amount of small 
scale power to produce the needed amount of bound ob- 
jects at the appropriate (high) redshifts Q. The issue 
is particularly important because the WMAP team data 
analysis seem to favor a running spectral index, when no 
link between reionization redshift and value of the run- 
ning spectral index is assumed. 

This difficulty may be overcome if the recombina- 
tion/reionization history is altered by some non-standard 
process, like black hole evaporation or the decay of heavy 
particles, and theorists have started to develop specific 
models to this aim |^. Bean et al. Q analyzed the im- 
pact of the recombination history implied by the mod- 
els. As for decaying particles, their analysis is limited to 
fairly long-lived specific Supermassive Dark Matter par- 
ticles. In a recent paper, Hansen and Haiman (here- 
after HH) have proposed a model in which the products 
of heavy sterile neutrino decay would significantly reion- 
ize the Universe at redshift z ~ 20. At difference with 
the previous models, these neutrinos may have a shorter 
lifetime and an inferred abundance linked to a different 
physics. 

Given the apparent variety of decaying particle mod- 



els, we take here a general approach and ask whether 
current data already put constraints on the reionization 
history implied by the decay and, in turn, on the parti- 
cle physics model involved. We revisit the reionization 
history implied by decaying particles and compare the 
implied CMB and large scale structure spectra with cur- 
rent data. We initially perform our calculations in the 
HH model, which we then generalize to the case of un- 
specified decay time and particle abundance. We work 
in a flat Universe with = 0.3, JIa = 0.7, Vlf, = 0.05 
and Ho = 70 Km/s Mpc^i. 

2) The decaying particle model for reionization 
We consider here particles with a decay time tdec of the 
order of 10^^ s (and define T15 = tdec/W^^s). The typical 
decay redshift for such particles would be Zdec < 20. For 
clarity, we refer here to the model presented by HH, in 
which a massive neutrino with mass of a few hundred 
MeV decays into an electron plus a pion. The main mech- 
anism of reionization in this case is the following: the 
relativistic electrons produced by the decay process in- 
verse Compton scatter the cosmic microwave background 
(CMB) photons, which then rcionize the hydrogen atoms. 
Each electron would be able to cause the reionization 
of a number of photons, according to its energy Ee = 
(M — l)m7r/2 where M = mx/m-n > 1. Typical electron 
energies are of the order of 1-100 MeV; each electron is 
therefore expected to result in about 10^ ionizations, with 
an efficiency of about 1/3. The sterile neutrino comoving 
number density is (HH): — 5.5 x 10~'' ^^^j.j^^2_i) 
which would imply an abundance at the present time of: 
ilx — 7 X J M'^ , irrelevant for determining the 

expansion at any time |l3l |. 

The decay process produce an additional ionization 
source which should be added to the standard rate in 
the ionization equation: 

where F^; = dn^/dt in comoving coordinates, Ih = 
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FIG. 1: Different rates are plotted as a function of redshift. 
Witliin the "standard" decaying particle model {nix ~ 200 
MeV, ri5 = 4)), the solid line is the recombination rate in 
the decaying particle model, the dotted is the ionization rate 
due to the decaying particles and the straight line is the ex- 
pansion rate H(z). The higher ionization produced by the 
decay particle implies a higher recombination rate at redshift 
2 ~ 10 — 600 with respect to the standard reionization one 
(dashed line). 



13.6 eV and the efficiency e is approximately 1/3, leading 
to an energy per reionization Ai? ~ 50 eV. 

Because only a fraction of the electron energy is used 
to produce ionization, the remaining part is assumed to 
heat the intergalactic medium. The additional term in 
the equation for the matter temperature Tm therefore 
reads: 



(dTm/dt)^^ 



1 



Ih + 3/2fcBT„ 
(Ih/c) 



(2) 



We modified the RECFAST code 8] by adding the 
terms in eq. ^ and eq. Q • Fig. ^ shows the reioniza- 
tion/recombination rates for this ionization prescription. 
We have chosen as reference parameters nix = 200 MeV 
and ri5 = 4, "which would imply that most particles have 
decayed by z ~ 20. The presence of the extra- ionization 
source boosts the recombination rate at redshift z < 800 
(fig. This means that most of the atoms reionized 
by the decay process would quickly recombine. As a re- 
sult, the maximum ionization fraction produced with this 
reionization prescription never exceeds 0.03 (fig. in 
contrast with what claimed in HH who predict a 50% ion- 
ization fraction a; at z = 20. Within the sterile neutrino 
model, it is impossible to reach such level of reionization 
just changing the parameters. A T15 increase only shifts 



FIG. 2: Ionization histories with and without a decay model. 
The dashed line correspond to standard recombination with 
no reionization. The solid line refers to the standard decay 
model adopted here {rrix = 200 MeV, ns = 4); note that the 
maximum reionization implied is 2; ~ 0.03. The dot-dashed 
lines correspond to a decay model with the same particle mass 
but ri5 = 8 — 12, while the dotted lines correspond to the same 
ri5 = 4 and rrix = 350 — 500 MeV. The three dot-dashed line 
is obtained by the reference model by increasing the number 
density of the decaying particle by a factor 300. 



the maximum value of x toward lower redshifts, while al- 
tering the mass of the particle cause the neutrino number 
density and electron energy to have compensating effects 
on the ionization fraction. 

In order to obtain 50% of reionization at z = 20 as a 
consequence of the decay process, the number density of 
the neutrinos must be artificially increased by a factor of 
order 300 (see fig. 13). This is not feasible in the sterile 
neutrino model, but it may be conceivable that other par- 
ticle physics candidates would produce such abundance. 
It is therefore sensible to ask which is the maximum level 
of reionization that can be produced by decaying particles 
without violating the actual CMB constraints. If these 
particles produce yields which then cause the Universe to 
reionize, the reionization process is likely to be similar to 
the one described above. In order to investigate to what 
extent such reionization process is allowed, in the follow- 
ing we take as a reference model the rrix — 200 MeV 
and ri5 = 4 sterile neutrino case to compute a reference 
abundance and decay yield energy. We then analyze the 
consequences of an hypothetical particle with abundance 
fx times the one of the reference case {ilx = 2.34 x 10^^). 
We consider decay times in the range 1 < ris < 20. 

3) Implications on the CMB power spectrum In this 
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FIG. 3: Temperature power spectra of the CMB. The solid 
line is the model with no reionization. The dashed is a model 
with sudden reionization at 2: ~ 17 (implying the WMAP best 
fit value r — 0.17). The three dot-dashed is the standard 
neutrino decaying model from HH, and the dot-dashed line 
is the decaying model where the abundance of the particles 
have been increased by a factor 300 in order to have a:: ~ 0.5 
at z = 20. 



section we examine the impact of the modified reion- 
ization history on the CMB power spectrum. It is well 
known that reionization has different effects on the tem- 
perature (XT) power spectrum and on the temperature 
polarization cross-correlation (TE). In both cases, an in- 
creased optical depth r causes a damping of the spectra 
that is progressively more pronounced for higher Vs. The 
TE spectrum, however, is boosted at low Vs if significant 
reionization occurs at low redshift. The earlier reioniza- 
tion occurs, the greater the I's where the effect appears 
in the TE spectrum. While the overall optical depth is 
to a large extent degenerate with other parameters, the 
low~rs bump in the spectra is a quite unique signature 
of reionization, and is precisely what has been observed 
by WMAP. Assuming instantaneous and complete reion- 
ization occurring at some low redshift Zre, the WMAP 
team has found r ~ 0.17, corresponding to z^e — 17 0. 

In the neutrino decay reference model, the decay pro- 
cess implies an alteration of the reionization history al- 
ready at high redshift (z < 800). Despite the fact that 
the ionization fraction x never reaches 0.1, it remains sig- 
nificantly higher than in the standard recombination case 
for an extended period of time, during which the number 
density of baryons is also high. This causes the optical 
depth to z ~ 1000 to be high even if the ionization frac- 
tion never reaches high values. For the neutrino decay 



FIG. 4: Temperature-Polarization correlation for the models 
of fig. |3] The lower dashed curve is a model with sudden 
reionization ai z — 6, which falls short in fitting the low Ts 
power spectrum. 

model with irix = 200 MeV and T15 = 4, we find that 
the optical depth in the range 20 < z < 800 is t = 0.11. 
These models therefore provide a way of fitting the TT 
power spectrum at high I's without significantly alter- 
ing the current accepted ranges for other parameters (see 
fig. Note, however, that a significant level of reion- 
ization at z ~ 10 — 15 is needed in order to match the 
WMAP low-Z result of the TE power spectrum, and sud- 
den reionization at z = 6 would fall too short in fitting 
the low-/'s TE spectrum even if the adequate total opti- 
cal depth was mainly provided by the decaying particles, 
(fig- El- For this reasons, we still have to invoke an early 
star reionization. 

Can decaying particles produce a reionization fraction 
of 0.5 at z ~ 20 ? As shown in fig.El in such a reionization 
scenario the ionization fraction would always be a; > 0.02. 
The total optical depth would then be extremely high, 
causing an excessive damping of the TT and TE power 
spectra at high Vs. Moreover, the reionization signature 
in the TE power spectrum would appear at too high I 
values {I ~ 30) contradicting the data. We conclude that 
it is not possible to reionize the Universe to high levels 
at z ~ 15 with a process that is powered by dark matter 
particle decay. 

4) Cosmological and astrophysical constraints We 
want to determine to what extent current cosmological 
data allow a particle decay process with a decay time 
shorter than or equal to the age of the Universe to 
trigger reionization. We performed a multi-parameter 
fit to the most recent radiation and matter power 
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FIG. 5: Marginalized likelihoods for the decay model of reion- 
ization parameters. 



spectra data. We allowed for particle decay reionization 
with timescales 1 < ri5 < 140 and arbitrary abundance 
(parametrized by fx)- In addition, we allowed for 
instantaneous reionization at a redshift Zre- We kept 
the spectral index Us and the amplitude of the matter 
power spectrum As as free parameters. The results 
were obtained with a modified version of the CosmoMC 
package [l3| which uses the WMAP likehhood [HIi^. 

The results are shown in fig. |S1 The abundance of 
decaying particles is limited by < 51 at 95 % C.L., 
and there is no evidence in the data that such decay- 
induced reionization models fit the spectra better than 
the standard reionization scenario. Longer lived particles 
are preferred to short lived ones. The data naturally con- 
strain the reionization rate in the early Universe, which 
is proportional to /x/tiq. Its 95 % C.L. upper limit is 
fx/Ti5 < 0.48. The reason why long-lived particles are 
favored resides precisely in the fact that there is a broad 
range of possible abundances for which the reionization 
history is little affected. On the contrary, short-lived 
particles imply a higher reioniation rate, therefore con- 
straining fx more. 

The marginalized likelihoods show a peak at a slightly 
lower reionization redshift (z^e — 13) and a smaller opti- 
cal depth with respect to the standard reionization case 
investigated by the WMAP team, but the likelihoods are 
still wide. The spectral index is not significantly affected. 

As for other astrophysical constraints, it has been ar- 
gued (HH) that decay-induced reionization may produce 
a specific distortion in the cosmic microwave background. 



We computed the evolution of matter temperature with 
the extra-term in eq. [5] In the decay particle scenario 
the additional ionization causes the matter temperature 
to follow the radiation one for a longer time. The de- 
rived y parameter, which depends on the difference be- 
tween matter and radiation temperatures in the past, is 
therefore very small. For the reference neutrino model, 
it is \y\ < 3 X 10^^°, and increasing the decay particle 
abundance by a factor 300 we have < 4 x 10^^"; still 
too small to be detected. 

5) Conclusions In this paper we analyzed to which ex- 
tent decaying particles may contribute to produce the 
high ionization rate observed at high redshift by the 
WMAP experiment. 

We showed that if decay yields had to produce a sig- 
nificant reionization at redshifts z ~ 20, they would alter 
the reionization history already at redshift z ~ 800, sig- 
nificantly modifying the reionization/recombination rate 
and the derived ionization fraction at all redshifts. The 
implied optical depth would be too high to match the 
one derived from recent CMB observations. Therefore 
the decay mechanism cannot be invoked to produce a 
significant ionization fraction at high redshift. 

As a consequence, a standard reionization mechanism 
(e.g. starlight) acting at fairly high redshifts {zre — 13) 
must still be invoked in order to fit the polarization data. 
These findings suggest that having some reionization pro- 
duced by decaying particles wouldn't ease the constraints 
put on small scale matter power spectrum by the require- 
ment of having enough collapsed structure at high red- 
shifts to produce small scale objects. Decaying particles 
do no alleviate the challenge that current data impose to 
models of structure formation like WDM or CDM with 
a running spectral index. 

We showed that the current CMB and large scale struc- 
ture data constrain the abundance of the decaying dark 
matter particles to < 1.3 x 10-'^{l + zf at 95 % C.L., 
when the decay time is in the range 1 — 140 x 10^^ s. 
Within this range, the data tend to prefer long-lived par- 
ticles. This result translate in a limit on the fraction of 
electrons reionized per unit time in the early Universe: 
2 X 10^^® s^^. There is no evidence that this altered 
reionization prescription produce a better fit to the data 
than the instantaneous reionization one. We evaluated 
the \y\ distortion that would be expected from this reion- 
ization mechanism, and found it to be below 10~^, there- 
fore too small to be observed. 
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